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We consider squeezing of one component of the collective spin vector of an atomic ensemble inside 
an optical cavity. The atoms interact with a cavity mode, and the squeezing is obtained by probing 
the state of the light field that is transmitted through the cavity. Starting from the stochastic 
master equation, we derive the time evolution of the state of the atoms and the cavity field, and we 
compute expectation values and variances of the atomic spin components and the quadratures of 
the cavity mode. The performance of the setup is compared to spin squeezing of atoms by probing 
of a light field transmitted only once through the sample. 
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I. INTRODUCTION 

Interactions between light and matter have several ap- 
plications within quantum information processing. When 
a light field interacts with a collection of atoms, the light 
and the atoms become entangled, and the state of the 
total system can no longer be written as a direct prod- 
uct of quantum states of the individual systems. As a 
consequence, if the light field is subsequently subjected 
to measurements, the state of the atoms will also be af- 
fected. This has, for instance, been utilized to entangle 
two atomic ensembles [E 0| and to teleport the state of 
a light field onto atoms [J] . It has also been suggested to 
generate various squeezed and entangled states of light 
and matter by sending a light field twice [4] or multiple 
times @ through the same atomic ensemble from differ- 
ent directions. 

The generation of entanglement between light and 
atoms may also be utilized to perform a quantum non- 
demolition measurement of one of the components of the 
collective spin vector of an atomic ensemble [1, H, H B 
B HO) EH) Ell- The measurements can reduce the un- 
certainty in the measured observable below the uncer- 
tainty of a coherent spin state, resulting in a squeezed 
spin state. Apart from the fundamental interest in gen- 
erating squeezed states, spin squeezing can improve the 
precision of measurements of, for instance, weak mag- 
netic fields [3, [13] . The strength of the interaction be- 
tween light and atoms is normally weak but can be en- 
hanced by placing the atoms inside an optical cavity as 
depicted in Fig.[Tl and in the present paper we investigate 
the performance of spin squeezing in a cavity compared 
to spin squeezing in free space. 

In spin squeezing experiments the usual initial state of 
the atoms is a coherent spin state, where all the atomic 
spins are oriented in the same direction, which we shall 
take as the x-direction. If the number of atoms is large, 
the a;-component of the collective atomic spin J = j^, 
where ji is the total spin of the ith atom, may be treated 
as a classical quantity ~ {Jx), and the commutator 
between the scaled spin components Xat = Jyl{f^{Jx)Y^'^ 
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FIG. 1: Experimental setup for probing of the ^-component 
of the collective spin of an atomic ensemble with electromag- 
netic radiation. A continuous laser beam linearly polarized 
in the a;-direction and with photon flux $ enters the cav- 
ity from the left. The interaction between the light and the 
atoms, described by the Hamiltonian H, rotates the polariza- 
tion vector of the light field by an amount, which depends on 
the z-component of the atomic spin. The angle of rotation can 
be measured by performing a detection on the light leaking 
out of the cavity (ki and K2 denote cavity decay rates). The 
half wave plate transforms the field operators ax and hy for x- 
and y-polarized light into (fla; -f ay)/%/2 and (ax — ay^l\f2, and 
these polarization components are subsequently separated by 
the polarizing beam splitter (PBS). The measurement out- 
come is the difference in photo current between the two photo 
detectors. 



and Pat = Jz/ {fi'{Jx)Y^'^ turns into the canonical com- 
mutator [iat,Pat] = i- In this approximation the ini- 
tial coherent spin state is a Gaussian state, and since 
the interaction Hamiltonian and the measurements trans- 
form Gaussian states into Gaussian states as long as 
Jx can be treated classically, a very efficient Gaussian 
formalism, which provides several analytical results for 
both pulsed and continuous wave fields, is applicable, as 
demonstrated for free fields in Refs. The Gaus- 

sian description is easily generalized to take an optical 
cavity into account [15], but although we shall be mainly 
concerned with the limit of a large number of atoms be- 
low, we also demonstrate that analytical results can be 
obtained even without assuming the Gaussian approxi- 
mation for the collective atomic spin. 
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The paper is structured as follows. In Sec.|TT]we apply 
the stochastic master equation for the setup in Fig. [1] to 
derive expressions for the time evolution of the state of 
the atoms and the cavity field, and in Sec. IIIII we evalu- 
ate the variances and mean values of the collective atomic 
spin operators and the quadratures of the cavity field as 
a function of time in the limit of a large number of atoms. 
Effects of losses due to spontaneous decay is considered 
in Sec. IIVI and the results are compared to those ob- 
tained for squeezing in free space. Section |V] concludes 
the paper. 

II. ATOMS INTERACTING WITH 
OFF-RESONANT LIGHT IN A CAVITY 

We consider atoms with a spin 1/2 ground state \g^) 
and a spin 1/2 excited state |e;p) interacting with a 
strong, off-resonant cavity field, which is initially lin- 
early polarized in the a;-direction. Decomposing the light 
into right and left circularly polarized cavity modes with 
field annihilation operators a+ — {—ax + iay)/^/2 and 
a_ ~ {cix -¥ idy) / ^/2, respectively, the Hamiltonian takes 
the form 

Wat 

H = hg^{a+\e+^^){g-^^\ + a_|e_,j)(g+,j| -h h.c.) 

i=l 

Wat 

-nA^(|e+,.)(e+,.| + |e_..)(e_,.|) (1) 

i=l 

in a frame rotating with the frequency of the light field. 
The summation runs over the TVat atoms, hg = —dEo, d 
is the atomic dipole moment, Eq = y^huj/{Veo), w is the 
angular frequency of the light field, V is the mode vol- 
ume, eo is the vacuum permittivity, and A = w — Wat is 
the detuning between the light field and the atomic tran- 
sition. For sufficiently large detuning g/A <C 1, the ex- 
cited states will not become significantly populated and 
can be adiabatically eliminated, which leads to the effec- 
tive Hamiltonian 

Wat 

i=l 

(2) 

In this and the next section, we neglect loss of photons 
and atomic coherence due to spontaneous emission from 
the excited states, but we return to an analysis of the 
role of spontaneous atomic emission of light in Sec. IIVI 
Since Y.^ \9^,^){9tA = iVat/2 ± J./?i, 

H = — — —(Kax + a'yay) - i— «ay - ala^) y. 

(3) 

The first term gives rise to a common phase shift of the 
X- and j/-polarized light and can be compensated by in- 
troducing an additional phase shift in the cavity, for in- 
stance by adjusting the length of the cavity. We thus 



ignore this term in the following. For the setup depicted 
in Fig. [2 it is desirable to have a large number of pho- 
tons in the x-polarized mode, since this increases the 
strength of the light-atom interaction, and since, in the 
polarization rotation measurement, the x-polarized field 
acts the same way as a local oscillator in balanced ho- 
modyne detection. Approximating d^ by its expectation 
value {dx(t)), the infinitesimal time evolution operator 
corresponding to the second term in Q takes the form 

U = exp ((5VA) {{dx{t))dl - {dUt))dy) {J,lh)dt^, and 

comparing this to the displacement operator D{5) — 
exp(5aj — (5*aj,), we observe that the interaction displaces 
the y-polarized field amplitude by an amount, which is 
proportional to the z-component of the atomic spin. De- 
tecting the quadrature of the j/-polarized field in the di- 
rection of the displacement thus constitutes an indirect 
measurement of Jz as stated above. 

Assuming a high finesse cavity kt ^ 1, where k is the 
total cavity decay rate and r is the round trip time of 
light in the cavity, and an only infinitesimal change of 
the atomic quantum state on the time scale of r, we de- 
duce from the detailed derivation in [l^ that the density 
operator p{t), describing the state of the atoms and the 
X- and y-polarized cavity fields, satisfies the linearized 
stochastic master equation 

dp{t) = - -^ [H,p[t)\ dt+^T]dK2 {dyp{t)dt + p{t)dl) dtjs 
+ ^ i~'^l"'xP{t) - p{t)dl.dx + 2dxp{t)d\) dt 

+ y^i3{t)[dl, p{t)]dt - ^p*{t)[dx,pit)]dt, (4) 

where k = ki + K2 + n^, ki (^2) is the cavity decay rate 
due to the lower left (right) cavity mirror in Fig. [l] is 
the cavity decay rate due to additional losses, /3(i) is the 
amplitude of the incoming probe beam, rj^ is the detec- 
tor efhciency, and dy^ is a stochastic variable representing 
the measured difference in photo current at time t (see 
[1^ for details). The first term in ^ is the Hamilto- 
nian evolution due to the interaction between the atoms 
and the cavity modes, the second term represents the 
knowledge obtained from the continuous measurement, 
the third and fourth terms take cavity decay into account, 
and the fifth and sixth terms appear due to the presence 
of the input beam. The derivation in assumes that 
y/Ki\P{t)\T is small, but, for a classical x-polarized mode, 
it is sufhcient to assume that ^^^^{t) varies slowly within 
times of order r (and if we are not interested in features 
of the solution occurring on a time scale r or faster, we 
may even allow y/l^(3{t) to change abruptly). In fact, if 
the x-polarized light is used as local oscillator as in Fig. 
[U it is required that 4KiK2|/?(i)pT/K^ » 1, since the 
local oscillator is assumed to be strong. We note that 
the stochastic term in does not conserve the trace 
of the density operator, which should hence be normal- 
ized explicitly. The probability to obtain the normalized 
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state p{t + dt) /Tr{p{t + dt)) at time t + di, given that the 
normahzed state at time t was p(t), is Tr(p(i + dt)) mul- 
tiphed by the probabihty to obtain the required value of 
dj/s, assuming that dys is a Gaussian distributed stochas- 
tic variable with zero mean value and variance dt. If it is 
the reflected light and not the transmitted light, which is 
subjected to measurement, and if the lower right cavity 
mirror in Fig. [1] is perfectly reflecting, we note that 'qdK2 
should be replaced by jy^Ki and k = ki + K2 + k^, should 
be replaced by k = ki + in Eq. (|3]). 

Since the Hamiltonian ([3]) commutes with , we can 
restrict ourselves to the basis consisting of the states with 
total spin quantum number J — -/Vat/2 if the initial state 
is a coherent spin state. We thus consider the states 
J,\n) = hn\n) with n ^ -N,t/2, -N,J2 + 1, . . . , 7Vat/2, 
and write the density matrix as 



p{t) = ^^p„™|?i)(m| 



(5) 



where pnm are operators on the space of the x- and y- 
polarized cavity field modes. This leads to the (iVat + 1)^ 
independent equations 



dl&x) Pnmdt 

^ • '-r nmi^^dy dy(ix)dt 

nm nm 
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-mp, 



2 



^y^y Pnm Pnmdydy 2,CLy Pnmdy^ dt 

~^ y djj-dx Pnm Pnnid\.Ojx 4^ 2(ix Pyimd^x) ^t 
+ y/K^l3it)[dl, p„rn]dt - y^P* (t)[ax, Pnm]dt (6) 

with solution 

ynm = Cnm{t)\-in{t))x{lm{t)\ ® Mt)) it)\, (7) 

where |7„(t)) and |a„(t)) are coherent states satisfying 



dt 



hi , , 

-^ln{t) 



n^anit) + y^Pit) (8) 



and 



da„{t) K g"^ 

= -2^n{t)+n—jn[t). 



dt 



For a classical input field the term in 
an{t) is negligible, and, assuming f3{t) = f3*{t) 
and a„(0) = 7ri(0) 0, we obtain 



(9) 

proportional to 



{dx{t)) = 7„(i) = 



,-^{*-t')/^P[t')dt' (10) 



and 



n{t) = n^-ir / e-''^'-''^'\dx(t'))dt' = na{t), (11) 



where a(t) is real and independent of n. Under these 
conditions the coefficients C„m in satisfy 

-p7^ = y^m^in + m)ait)dys--{n-m)^a{t)^dt (12) 

with solution 



C^m{t) = C„™(0) exp ( - 2 (" - aitlydtl 



+ y/Vdi^in + m) a{t')dy' 



-']^(^n + m)^J a{t')^dt'y (13) 

If required, the analysis is easily generalized to include 
all simultaneous eigenstates of and , since it is only 
needed to include more terms in ([5]) and to introduce 
additional labels to distinguish the different states. 

We finally note that if $(i) is zero for < < and as- 
sumes the constant value $ for i > and if the light- 
atom coupling is sufhciently weak to ensure that the 
change in the state of the atoms during the transient 
is negligible, we may approximate {dx{t)) and a{t) by 
their respective steady state values {dx) = 2\/ki$/k and 
a = 2g'^{ax) / (kA) for t > 0, which makes the integrals in 
trivial to evaluate. In that case the state at time t de- 
pends on the measurement result through the integrated 
signal Ys = /g dys only, and the probability density to 
measure a given value of Yg is [lB| 



P 



C„„(0) 



■ exp 



{Ys - 2^r]dK2natf 
2t 



(14) 



The measurement leads to a narrowing of the distribution 
Cnn{t) / Y^.^Cmm{t) ovcr the possiblc eigenstates of J^, 
but the expectation value of depends on Yg , and if we 
average over all possible measurement outcomes, we find 
that C„„(i) = C„„(0). 



III. EXPECTATION VALUES OF ATOMIC SPIN 
OPERATORS FOR LARGE A^at 

Having obtained the state of the atoms and the y- 
polarized cavity field as a function of time, we can 
now evaluate expectation values and variances of the 
atomic spin operators and the field quadrature opera- 
tors Xph = [dy + dl)/^/2 and pph = -i{ay - at)/\/2. 
We assume below that the initial state is a coherent spin 
state pointing in the x-direction and that the number 
of atoms is large iVat ^ 1, since this is a typical ex- 
perimental condition, and since it allows us to simplify 
the obtained expressions considerably. In order to stay 
within the parameter regime where the Gaussian approx- 
imation, discussed in the Introduction, is valid, it is also 
required that the total measurement time is short com- 
pared to the time it takes to gain sufficient information 
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to project the state of the atoms onto a single eigenstate 
of Jz- For the steady state case it follows from Eq. 
that the relevant time scale is determined by the condi- 
tion AridK,20?t ~ 1, and we thus assume in the following 
that 4?7(i/t2 j^a{t'Ydt' is small, i.e., comparable to the 
size of N~^^, while 2^ridH2 a{t')dy'g is assumed to be 

— 1/2 

comparable to 

First we would like to determine whether the atomic 
spin is indeed squeezed, and we thus trace out the cavity 
field and compute the variance of 



For a coherent spin state pointing in the x-direction 



(15) 



Cnm (0) 



1 



iVat! 



2^- V iN.t/2 + ny.iN,t/2-ny. 



y (A^at/2 + m)!(iVat/2-m)!^ 
we may apply the approximation 



Cnm(O) 



■ exp 



(16) 



(17) 



TTiVat V ^at 

and it follows from ((13]), (US]), and ^7} that 
YaT(j,/h) 1 / ft N-i 

-^tj^--2[^+^^^^^-^l-(^rdt') . (18) 

Remarkably, this result does not depend on the measure- 
ment re adout and is thus deterministic. The variance of 
Jz/h/ is seen to be decreasing and smaller than 
1/2 at all times t > if rjd > 0. For the special case where 
<i>(t) is zero for < < and assumes the constant value $ 
for t > 0, we have {axit)) = y^4ki$/k2(1 - exp(-Kt/2)), 
and 



Var 



^ ^^ ^AT 4ki$ 4/ r]dK2 - 



where 



t = t- — 



11 2Kt + ?> 



-Kt/2 



n'^f + 6Kt + 10 



2k k 4:K 
This is to be compared to the expression 



(19) 



(20) 



/VarGVfi)^ 



1 / 

\ / sp 



for single-pass squeezing jlll |. Apart from what ef- 
fectively amounts to a small reduction of the probing 
time, appearing because it takes a short while to build 



up the cavity field, the effect of the cavity is to in- 
crease the coefficient multiplying t by a factor Q — 
16kiK2/('*^''"^). In the single-pass case each segment 
of temporal width r of the probe beam interacts only 
once with the atoms, and (cix) ~ V $t for all times 
t > 0. The interaction thus transforms the y-polarized 
mode from the vacuum state |0) into the coherent state 
;7|0) = \n{g'^T/ A)V^), where, for simplicity, we have 
assumed that the atoms are in the Jz eigenstate \n). The 
number of y-polarized photons observed per unit time is 
thus n^{g'^T^ / A^)^'rid. If the cavity is included, on the 
other hand, the number of y-polarized photons observed 
per unit time is the product of the number of y-polarized 
photons in the cavity |anp, the rate K2 with which the 
photons leave the cavity through the cavity output mir- 
ror in Fig. [U and the detector efficiency r]d, and the re- 
sult n^(45^/(K^A^))(4Ki<I>/K^)r7dK2 is larger than in the 
single-pass case by precisely the factor Q. To under- 
stand this increase in the number of detected y-polarized 
photons, Q may be divided into the three factors K2/K, 
4Ki/(K^r), and 4/(Kr), where the first appears because 
the effective detector efficiency is ridK2/K for squeezing 
in a cavity and rjd for single-pass squeezing, the second 
factor appears due to the increase in the number of pho- 
tons in the a;-polarized mode, as can be seen from the 
increase in production rate of y-polarized photons, when 
the flux of x-polarized photons in the case of single-pass 
squeezing is increased from $ to 4ki$/(k^t), and the 
third factor appears because photons are present in the 
y-polarized mode in the cavity, as can be seen by compar- 
ing the number of produced y-polarized photons when U 
acts on |q;„) and when U acts on |0). 



Applying Cn±in^iiO) 
we also find 



C„„(0)(l-2/7V^t + O(iV,r)), 



Var 



iVat/2 



i ( 1-l-iVatK 



a{t'fdt' + iVata(t)' 



(22) 

Since k is larger than or equal to ridi^2, the product of 
and ([^ is larger than or equal to 1/4 as required 
by the Heisenberg uncertainty relation. Equality is only 
obtained for a{t) = and k — r/dK2, where the first equa- 
tion is satisfied if the y-polarized cavity mode is in the 
vacuum state at the final time t, and the second equa- 
tion is satisfied if all photons that leave the cavity are 
detected. 

It follows from 



Var(xph) = - - 



and 



N^tajt)^ + N,tVd'i2 Jo ajt' fdt' 
1 + A^at%K2 /o a{t'Ydt' 



(23) 



Var(pph) = 1/2 (24) 

that the cavity field is not squeezed, but, for time in- 
dependent a{t), the uncertainty in iph decreases with 
probing time. The Heisenberg limit is only achieved ex- 
actly at times, where the cavity field is in the vacuum 
state. 
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The expectation values 

{Jy/h) = 



2/iVat + 2r/rfK2 /o a{t'Ydt' 
V2ait){.h/h) 




(25) 
(26) 

(27) 
(28) 



are either stochastic or zero, depending on whether the 
measurements supply information on the concerned oper- 
ator or not. Different mean values of Jz are thus obtained 
if the experiment is repeated. By applying feedback and 
rotating the collective spin, it is, however, possible to 
achieve absolute squeezing, where the same mean value 
of Jz is obtained in each run [l^l . 

An alternative approach to calculate expectation val- 
ues and variances of Xat — Jy/{^{Jx)Y^'^, Pat — 
Jz/{h{Jx)y^^ , Xph, and Pph is to assume from the start 
that the state of the atoms and the y-polarized cavity 
mode is approximately Gaussian at all times t satisfying 
477^^2 /p a(i')^(ii' <C 1. Gaussian states are efficiently de- 
scribed in terms of Wigner functions, and we thus trans- 
late the nonlinear stochastic master equation for the den- 
sity operator for the atoms and the y-polarized cavity 
mode derived in uM 



dp{t) = --[H,p{t)]dt 

+ iap{t) - Tr (ap(t)) p{t)) dWs 

+ {Pit)^^ - Tr {p{t)a)) p{t)) dWs 

+ {-a)ap{t) - p{t)a)d + 2dp{t)a)) dt, (29) 

where dWs is a Gaussian distributed stochastic variable 
with zero mean value and variance dt, into an equation 
involving the Wigner function W 



dW 



-g{t) Pat 



d 



' dx 



ph 



Pph 



d 

9Xa 



Wdt 



+ \/'^VdK2 ^a^ph - {xph) + WdWs 



1 



4 \dxl^^ dp^ 




Wdt, (30) 



where is a function of t and the quadrature variables 
a; at, Pat, Xph, and pph, and we have introduced the effec- 
tive light atom coupling strength 



2^2 /(J,) (a,(0) 
in terms of which the Hamiltonian reads H 



(31) 



hg{t)pa,tPph- For a Gaussian state 
1 



W 



7rVdet(F) 



exp{^iy^{y)fV-\y-{y))), (32) 



where y — {xat,Pat,Xph,Pph)'^ is a column vector of 
quadrature variables, y = {xat,Pat,Xph,Pph)'^ is a col- 
umn vector of the corresponding quadrature operators, 

andv = {{y-{y))iy-{y)r) + {iy 

the covariance matrix. Inserting 
that 



T\T 



IS 



(ym-m 

into (PU)) . we find 



where 



dV_ 
'dt 



G 



G - DV -VE ~ VFV, 



D = E^ = 



- 





























K - 77K2 





. 








K 























(33) 



(34) 



-g{t) k/2 - rjdK2 














k/2 



F = 



" 








0" 




















ridK2 





. 












(35) 



(36) 



and 



d{y) 



- 





m 











m 


-k/2 





. 





-k/2 



{y)dt 



(V-I) 



dWs. (37) 



Equation (j33|) is a so-called matrix Ricatti equation, and 
if V is decomposed according to ^ = AIK~^, it can be 
rewritten as the linear set of equations M = —DM + GK 
and K = FM+EK. Solving these equations analytically 
for a time independent g, we find expressions, which are 
in accordance with the above results. Equation p3p can 
also be derived following the covariance matrix approach 
outlined in Ref. [Ill for single-pass interaction. To do 
so, the light beams are divided into segments of duration 
T, where each segment constitutes a classical z-polarized 
field mode and a quantum mechanical y-polarized field 
mode, and the state of the atomic spin and the quantum 
mechanical field modes is assumed to be Gaussian. The 
time evolution of the covariance matrix is then obtained 
by realizing that an interaction between the atoms and 
the field modes, a beam splitter operation, and a ho- 
modyne detection of a field mode all amount to certain 
transformations of the covariance matrix. 
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IV. INCLUSION OF LOSS DUE TO 
SPONTANEOUS DECAY 

If the atoms are allowed to decay by spontaneous emis- 
sion, there will be a loss of atomic coherence as well as 
a decay of the mean spin, because the polarization of a 
spontaneously emitted photon, in principle, provides in- 
formation on the final state of the atom that emitted the 
photon. To include spontaneous emission in the analysis, 
we add decay terms to the master equation for interac- 
tion of the atoms with an a;-polarized and a y-polarized 
light mode 



dp 



r( -D[\g^,,){e+4]p^ 



'-D[\9^ 



+ ^i^[|5+,^)(e-,.|]p-Hii^[|5-,»>(e- 

where H is given by ([l]) and 

D[c\p = cpc' — {c^cp pc)c)/2. 



(38) 



(39) 



Adiabatic elimination of the exited atomic states leads 
to 



dp — 



4 



, r 

dt+-- 




'-t 4. 

+ 3a+l5-,»>(5-,j|p|5~,j)(5'-,»|a+ 
2 

+ 3a-l5-,»)(.9+,j|p|5+,»)(5-,»|a- 
+ \a+\g+,,){g^4p\g^^,){g+,,\d^^dt, (40) 

where, as before, g/A ^ 1 and we have omitted the term 
in the Hamiltonian giving rise to a common phase shift 
of the light modes. Finally, homodyne detection, cavity 
decay, and the input beam are taken into account by 
adding the terms 

,/m^{dp{t) - Tr (dp{t)) p{t)) dWs 

(P(i)a^ - Tr {p{t)d^) p{t)) dWs 
+^ {^°l°-yPi^) ^ P{i)°'l^y + 2ay/o(t)a],) dt 
+1 (-4axP(i) - p{t)ald^ + 2d:cp{t)dl.) dt 
+^(3{t)[dl,p{t)]dt- ^I3*{t)[d^,p{t)]dt (41) 



on the right hand side of pO)) . 

Equation ([40|) . (|4T|l can be solved numerically for a 
small number of atoms and a classical z-polarized mode, 
but here we aim at an approximate description, which is 
valid for the case, where the a;-polarized mode is classical, 
the initial atomic state is a coherent spin state pointing in 
the x-direction, Nat is sufficiently large to assume that 
Jx is classical, and t is small compared to the time it 
takes to project the atomic state onto an eigenstate of 
Jz due to measurements and small compared to the time 
it takes {Jx) to decay significantly. From the stochastic 
mater equation it follows that 



A ((fly - dl)Jy) 



I djJxjt)) ,,,, 9 



iUt)) 



2A2-i-r2/4 h 



(42) 



The ratio between the last and the first term is approx- 
imately {dx{t)){Jx{t)y^^T/{2A), which evaluates to 10^ 
for the parameters given in the caption of Fig. [2l and we 
thus skip the first term and obtain 



{Jx[t)) = — — exp 



(43) 



where we have defined the time dependent decay rate 77 (i) 
of the atomic spin as 



77(t) = {dx{t)y 

Similarly, for (a^(t)) we find 



g 



2 A^+r^/A' 



(44) 



d{ax{t)) 
dt 



K + e 



(45) 



where we have defined the photon absorption rate as 



2 A2 +r2/4' 



(46) 



We can now use the stochastic master equation to derive 
expressions for the time derivative of the first and sec- 
ond order moments of Jy, J^, iph, and Pph, and we find 
that, apart from third order moments appearing in the 
stochastic terms of the equations for the time derivative 
of the second order moments, these expressions contain 
only first and second order moments. Since the state 
of the atoms and the light field is nearly Gaussian under 
the above conditions, we approximate the third order mo- 
ments by a sum over products of first and second order 
moments to obtain a closed set of equations. We also ap- 
proximate V12, Vi3, V24, and V34 by zero, because these 
covariance matrix elements are zero if spontaneous emis- 
sion is neglected, and because the rest of the covariance 
matrix elements only couple to V12, V24, and V34 
through terms that are proportional to the small factor 
{T/2)g'^ /{A'^ + r2/4). Within these approximations we 
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0.4 0.6 
(/ ms 



FIG. 2: Uncertainty in pat as a function of time with atomic 
decay included (solid curves) and excluded (dashed curves). 
The upper and lower curves represent squeezing of the same 
atomic system in free space and in a cavity, respectively. Note 
the different time scales. The parameters are (see [IH): Nat = 
10^^ <E> = 10" s~^ for t > 0, A = 2 cm^ r = 3 • 10"^° s, 
A = 27r • 10^" Hz, A = 852 nm, T = 3.1 • 10^ s~\ d = 2.61 ■ 
10"^^ Cm, K = Ki = 27r ■ 3 ■ 10^ Hz (i.e., we observe the 
reflected light and assume kl =0), and rjd = 1- 



find that the time evolution of the covariance matrix is 
given by the Ricatti equation (|33p with 



G = 



{Mt}) 



ij{t) 

l-qit) 











K + e — 

K+e 



(47) 



D = E'^ = 



ri{t)/2 

T]{t)/6 

-g{t) {K + e)/2^ 





-m 




(K + e)/2 



(48) 



and F given by Eq. (|36p. Apart from a factor 1/3 in D22 
and £^22 and a factor 2/3 in G22, which appear as a direct 
consequence of the factors 1/3 and 2/3 in Eq. (|38p . this 
is exactly what is obtained by generalizing the Gaussian 
treatment of spontaneous decay in Refs. j^, [ll|, and [l^l 
to squeezing in a cavity. 

Integrating the Ricatti equation numerically, we obtain 
the lower curves in Fig. [2l where Apat = (Var(pat))"^''^- 
For the chosen parameters g{T <^ t r]{t ^ t)~^)t = 
2 • 10-3, KT = 6 • 10-3, and $r = 3 • lO'', and the re- 
quirements of a dilute atomic g high finesse cavity. 



and a strong local oscillator are satisfied. The values 
{X{0))/h = 5 • 10" » 1 and {a^{t > t)) = 4.6 • 10^ > 1 
justify the classical treatment of these quantities, and 
t = 1 ms satisfies t < (477dK2a(i > r)^)-^ = 2.7 • 10"* s 
and t ^ ri{t » r)^^ = 0.13 s. When atomic decay is 
included, the uncertainty in pat does not decrease indefi- 
nitely, but begins to rise at a certain point if the probing 
is continued. For the given example, the minimum value 
of the uncertainty is (Apat.)min = 0.0233. 
For single pass squeezing |ll| 



dVar(pat) 
dt 



4 2 

-27Vat$^%Var(pat)'e-''* 



A2 



-77Var(pat) 



(49) 



and the result of an integration of this equation is shown 
in Fig. [2] for comparison. The squeezing is seen to occur 
on a significantly slower time scale, and we note that 
Q = 5 • 10^ for the chosen parameters. The attained 
minimum value of the uncertainty (Apat)min — 0.121 is 
also significantly higher. This value is in accordance with 
the value 0.118 obtained from the approximate relation 



(Apa 



TyA^ 



1/4 



(50) 



derived in (we have included an additional factor 
of 2/3 to take the factors 1/3 and 2/3 in Eq. ([38]) into 
account). Since we found in Sec. IIIII that the main ef- 
fect of the cavity is to increase the squeezing rate by 
Q, and since it follows from ((44|) that 77 is a factor 
4ki/((k -|- e)2r) larger for squeezing in a cavity than for 
single-pass squeezing, we expect that (Apat)min is de- 
creased by a factor {{k + e)'^T/{4:K2)y^^ if the atoms are 
enclosed in a cavity. This leads to the predicted value 
(Apat)min = 0.0230 for Squeezing in a cavity, which is 
close to the value observed in Fig. [21 Since (Apat)miii is 
proportional to N^^^^'^, we could also regard the squeez- 
ing enhancement factor as a multiplicative factor on A^at i 
and this opens the way to use the cavity to achieve 
measurement induced squeezing of a smaller number of 
atoms. For A^at = 10^^ • (k + e)V/(4K2) = 1.4 • 10^ we 
thus find a minimum uncertainty of 0.121 after a probing 
time of 7 ms. We note that g{t), '7(0; and e are all un- 
changed if $, A'at, and A are scaled by a common factor, 
and we thus obtain the same result for 7 • 10^ atoms if 
$ = 5 • 10^^ s^^ and A = 10^^ m^. A further decrease in 
$ would, however, violate the assumption of a classical 
x-polarized field and the approximation below Eq. (|42|) . 



V. CONCLUSION 

We have considered squeezing of one component of the 
collective spin of an atomic ensemble achieved by per- 
forming homodyne measurements on light, which has 
interacted with the atoms, and we have found that 
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the squeezing rate can be increased by a factor Q = 
I6K1K2/ ('*^''"^) by placing the atoms inside an optical cav- 
ity. For ensembles containing a large number of atoms 
initially prepared in a coherent spin state, an efficient 
Gaussian formalism is applicable, from which we have 
derived equations for the time evolution of the covari- 
ance matrix describing the state of the atomic spin and 
the cavity field, but we have also demonstrated that an- 
alytical results for the state can be obtained even if the 
state of the atomic spin is not Gaussian. Despite the 
stochastic nature of the measurements, the variances of 
the components of the atomic spin and the quadratures 
of the light field evolve deterministically in the Gaussian 
approximation, and, in the lossless case, the variance of 



the squeezed atomic spin component is a monotonically 
decreasing function of time. According to the Heisenberg 
uncertainty relation the variance of the conjugate atomic 
spin variable has to increase, and the uncertainty prod- 
uct only attains the smallest allowed value if all photons, 
transferred to the mode with polarization orthogonal to 
the polarization of the probe beam due to the interaction 
with the atoms, have left the cavity and been detected 
at the considered time. Allowing the atoms to decay 
spontaneously, we find that the minimum variance of the 
squeezed spin component is obtained much faster and is 
approximately reduced by a factor ((k -|- e)^T/(4K2))^/^ 
compared to the single-pass setup. 
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